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Abstract 

The presence of ordered structures such as hehces in cohapsed states of polymer chains is stih 
an open question chahenging physics and biology. In this work, we present a potential model for 
polymer chains with monomers that are not strictly attractive, but that can make adhesive contacts 
with other monomers. We find that the chain develop helical order during the process of collapsing 
from an initially stretched conformation. It seems in this case that the adhesive contacts help the 
polymer chain to stay trapped in the helix state. 

PACS numbers: 36.20.Ey,82.37.-j,87.15.He 
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INTRODUCTION 



The three dimensional structure of a polymer chain will determine many of its properties. 
The most important examples of polymers are synthetic homopolymers, such as polyethylene 
an biopolymers, like proteins and DNA. In the case of synthetic polymers, the structure is 
the basis of a multitude of macroscopic products of everyday life[l|. On the other hand, 
biopolymers are the basic building blocks of cells and involved in the majority of biological 
processes such as chemical reaction catalysis, oxygen transport or muscles contract ions {2!. 
The ground states of these molecules are often crucial to their functions. 

In recent years, designing functional polymers at nanometers scales has gained increasing 
interest jsl. However, it is challenging to control structures that have the ability to undergo 
cooperative transitions between random conformations (globules) and ordered conforma- 
tions (a— helices and /3— sheets). There has been a spate of theoretical studies dealing with 
this issue, commonly achieved using Monte Carlo and molecular dynamics simulationsj4|. 
These works have offered rewarding insights for understanding the dynamics and the phase 
transition of polymers and there is still a big interest in investigating them until today ^. 

In a previous work, we have found that flexible homopolymers spontaneously develop 
helical order during the process of collapsing from an initially stretched conformation. We 
have also demonstrated that the helices are long lived transient states at low temperaturesjol, 
but with "borrowing" energy from their surrounding solvent particles, they overcome the 
energetic barrier and collapse into a stable globule. 

Here, we extend our study by adding to the polymer model adhesive contacts between non 
covalently linked monomers. These kind of contacts appears for example in the vital part 
of the inflammatory response, where the activated cells, due to adhesive contacts develop 
capabilities to transmigrate into the inflamed tissue 

The report is organized as follows. In Sec. II we describe the model, simulation method 
and the parameters used to quantify the ordered structures. In Sec. Ill, we present the 
results from the simulations. Finally, we conclude with a summary of our results and the 
appropriate acknowledgments. 
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MODEL AND METHODS 



In the simulation, we have used a simple bead-spring model. The polymer chain consists of 
identical monomers where each monomer position varies continuously in three dimensions. 
Dimensionless units are used during the simulation and are defined in terms of the bead size 
a and the Lennard- Jones energy e. There are three types of interactions: a harmonic bond 
length interaction between adjacent monomers 



with a = 100e/(T^ and Tq = 0.85a 

A Lennard- Jones potential between non adjacent monomers which is repulsive up to a 
distance x = a. 



In this potential, we have removed the attractive part but the excluded volume part is 
kept and adhesive interactions between any two monomers along the chain are possible. 
These interactions have three properties: 

Each monomer can only have one adhesive bond, i.e., if a monomer establish an adhesive 
contact with another monomer along the chain, establishing adhesive contacts with other 
monomers will be prohibited. This is important because it makes all the difference between 
the adhesive bond and general attractive interactions. To get collapse, the minimum energy 
of the bond should be negative. To study unraveling dynamics of the polymer chain, there 
should be an energy barrier between the adhesively bonded and non bonded state. 

n 

Assuming the validity of the Kramer's rate theory[8|, we have introduced an additional 
variable Xij that describes the state between two monomers i and j. 
If Xij = 1; there is an adhesive bond , if Xij = 0, there is none. 
The interaction potential for non-bonded monomers is then: 




(1) 
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with 
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FIG. 1: The non-bonded interaction potential 



A;((r/fr)2- 1.) -£ for r < 2.5(T 

Kdbond(r) = { (4) 

Vo otherwise, 

This potential is constructed such that the minimum value of Venter ~^ "^^^^ 
of bonded monomers (zero for non-bonded). Hence, an adhesive bond has the energy —e. 
K) = Knter(2.5(T) Is thc hclght of the energy barrier, thus the cutoff r^.^^ should be chosen 
large enough that Vq is positive. Bonds can only be established between monomers with 
distances less than r^.^ and may not be established between monomers that are neighbors 
along the chain (monomers n and n + 1), i.e., Xn,n+i = 0. 

In addition, there is a constraint "^jXij ^ 1? ^-S-) every monomer can have only one 
bond. 

In simulating the dynamics of this model. We have used a mixture of Monte Carlo for 
bond-breaking/bond-establishing and Langevin dynamics. One Monte Carlo step consists of 
enumerating all monomer pairs which are not bonded along the chain, than picking randomly 
one pair and if Xij = 1? break adhesive bond with probability exp(— /3(Vo — V{rij))). 

If Xij = ^i^d the bond is not forbidden by any constraint, establish adhesive bond with 
probability — /3Vo, than repeat this a number of times {e.g., Ntimes, once for each monomer). 
Langevin Dynamics step without hydrodynamics describes the motion of the monomers. 

Order parameter 

To characterize the helical structures, we use an order parameter that is commonly used 
to represent the net helical growth in chainlike molecules(9(|, as defined by 
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